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ABSTRACT 

Families  of  periodic  sequences  with  small  side-lobes  and  low  cross- 
talk (cross-correlations)  are  needed  in  many  digital  communication  systems. 
L.  R.  Welch  has  recently  established  a lower  bound  for  the  largest  of  the 
side-lobes  and  cross-correlation  coefficients  associated  with  a family  of 

M distinct  sequences  of  period  L.  In  particular,  for  two  sequences 

-k 

there  must  occur  a coefficient  greater  than  (2L)  . As  the  number  M of 

k 

sequences  increases  to  approximately  L , this  lower  bound  increases  to 

-k 

L . In  this  report  several  sets  of  periodic  sequences  utilizing  complex, 

ternary,  and  binary  coded  signal  bits  are  constructed  which  nearly  meet 

Welch's  bound  for  M approximately  equal  to  L . In  this  paper's  context, 

ternary  coding  refers  to  three  distinct  phase  signals,  while  complex  coding 

refers  to  phase  levels  relating  to  the  nC^  complex  roots  of  unity.  Host 

k 

of  these  sequences  contain  only  L non-zero  entries.  However,  one  type 
of  family  consists  of  M = 2m  binary  (±1)  sequences  of  period  L * 2™  - 1. 
The  maximum  coefficient  for  this  family  is  approximately  one-half  as  large 
as  the  corresponding  maximum  for  the  Gold  sequences  of  the  same  period. 

This  represents  a considerable  Improvement  over  the  Gold  codes  for  such 
codes  requiring  an  even  number  of  shift  register  stages  for  generation. 
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1.  INTRODUCTION 

The  need  for  sets  of  sequences  possessing  low  periodic  cross  correla- 
tions and  auto-correlation  side-lobes  frequently  arises  in  communications 
problems.  A recent  paper  of  Welch1  presents  a lower  bound  for  the  maxi- 
mum absolute  value  of  these  correlation  coefficients  as  a function  of  the 
period  L and  the  number  M of  sequences  in  the  set.  The  Gold  sequences,2 
which  have  period  L B 2m  - 1,  differ  from  the  bound  W(L,M)  of  Welch  by 
factors  of  approximately  2 and  2 for  m even  and  m odd,  respectively. 
Here  we  describe  sets  which  achieve  the  bound  W(L,M)  as  well  as  sets 
which  approach  it  as  L increases. 

The  sets  presented  in  Sections  2 and  3 are  formed  by  combining  partial 

difference  sets  in  cyclic  groups  with  Fourier  or  Hadamard  matrices.  A 

partial  difference  set  is  a slight  generalization  of  a planar  difference 

set,  which  gives  rise  to  a finite  cyclic  projective  plane. 3»l*  The  planar 

difference  sets  are  used  to  construct  optimal  and  near-optimal  sets  of 

sequences.  Additional  near-optimal  sequences  are  formed  using  certain 

of  the  relative  difference  sets  of  Flllott  and  Butson.5  With  only  one 

exception  the  optimal  sequences  have  complex  entries.  For  each  of  these 

2 

sets,  L is  approximately  equal  to  K , and  M equals  K or  K + 1,  where 


3 L.  R.  Welch.  "Lower  Bounds  on  the  Maximum  Cross  Correlation  of 
Signals,"  IEEE  Trans.  Inform.  Theor.,  May  1974,  pp.  397-399. 

2 R.  Gold.  "Optimal  Binary  Sequences  for. Spread  Spectrum  Multi- 
plexing," IEEE  Trans.  Inform.  Theor.,  October  1967,  pp.  619-621. 

3 L.  D.  Baumert.  Cyclic  Difference  Sets.  Berlin,  Heidelberg,  New 
York,  Springer-Vcrlag,  1971. 

• * M.  Hall,  Jr.  Combinatorial  Theory.  Waltham,  Mass.,  Blalsdell, 

1967. 

5 J.  H.  E.  Elliott  and  A.  T.  Butson.  "Relative  Difference  Sets," 
III.  J.  Math.,  Vol.  10  (1966),  pp.  517-531. 
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K Is  the  number  of  elements  In  the  underlying  partial  difference  set, 
and  L is  the  order  of  the  cyclic  group.  The  largest  correlation  coef- 
ficient has  magnitude  1/K,  and  each  sequence  has  exactly  K non-zero 
entries. 

The  sequences  described  in  Section  4 are  derived  from  binary,  cyclic 
error-correcting  codes  in  much  the  same  way  as  are  the  Gold  sequences. 

For  L = 2^m  - 1,  a certain  cyclic  error-correcting  code  of  length  L 
and  dimension  3m  yields  a near-optimal  set  of  2m  binary  sequences. 

Like  the  Gold  sequences,  these  have  no  zero  entries  and  are  easily  gen- 
erated by  shift-register  techniques.  Their  correlation  bounds  are  easily 
obtained  from  the  code  weight  distributions  given  by  Kasami,  Lin,  and 
Peterson.6 


2.  DEFINITIONS  AND  NOTATION 

12  M 

We  are  interested  in  sets  A = {a  , a , . . . , a } of  sequences  of 

period  L.  Each  aV  is  a complex  L-vector  •••  * ,)  of  norm  1; 

V v v v u 

that  is,  aQ  Iq  + ...  + aL  i i = 1,  where  a denotes  the  complex  con- 
jugate of  a.  The  correlation  coefficients  are  given  by 

cv»(T)  -Z  Vw  . 

i*0 

with  subscripts  reduced  modulo  L,  1 < V,  A < M,  0 < T <L-1.  A measure 

of  the  correlation  quality  of  the  set  A is  c^nv.(A) , the  maximum  over 

all  Icvx(t)I  wlth  V ^ A or  x ^ 0;  i.e.,  the  maximum  over  all  of  the 

auto-  and  cross-correlation  coefficients  except  the  "peaks”  c (0)  ■ 1. 

VV 

A principal  result  from  Welch1  is  that  for  any  such  family  A, 
c (A)  >W(L,M),  where 

mmX 


* T.  Kasami v S.  Lin,  and  W.  W.  Peterson.  "Some  Results  on  Cyclic 
Codes  Which  Are  Invariant  Under  the  Affine  Group  and  Their  Applications," 

Jnfom.  and  Contr. t Vol.  11  (1968),  pp.  475-496. 

4 
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wa;M>  ■ 

We  define  an  (L,K)  partial  difference  set  to  be  a set 

A - {d dR)  of  K distinct  elements  from  the  cyclic  additive 

group  of  integers  modulo  L,  such  that  for  every  non-zero  x in 
there  is  at  most  one  pair  df,  dg  from  A satisfying  d^  - dg  « x. 
From  this  condition  one  can  easily  show  that  L must  be  greater  than 
K2  - K.  Associated  with  the  partial  difference  set  A is  the  binary 
sequence  z(A)  = (z^,  ...  , z^  defined  by 


z 


i 


1 

if 

i G A 

0 

if 

i £ A 

The  requirement  that  d^  - dg  = x for  at  most  one  pair  from  A is 

equivalent  to  forcing  the  x1*1  auto-correlation  coefficient  of  z(A)  to 

be  0 or  1,  for  all  x =£  0.  If  all  the  side-lobes  are,  in  fact,  equal 

to  1,  then  A is  a planar  difference  set.  In  this  case  the  family  of 

translates  A + t,  t G Z , form  the  L lines  of  a cyclic  projective 

L 

plane  of  order  K - 1,  (see  Baumert3  and  Hall4).  (A  planar  difference 
set  is  a special  case  of  a cyclic  difference  set.3  For  the  general 
cyclic  difference  set  A,  the  side-lobes  for  the  associated  sequence 
*(A)  must  be  constant,  but  may  be  greater  than  1). 

Suppose  A is  an  (L,K)  partial  difference  set  and  E = (e^)  Is  an 

H X K complex  matrix  with  each  row  of  norm  I._  We  define  a set  A(A;E) 

of  H sequences  of  period  L as  follows.  Each  sequence  aV  is  zero 

except  for  the  entries  corresponding  to  elements  of  A.  The  K non-zero 
v th 

entries  in  a are  determined  by  the  v row  of  E.  More  precisely 
evr  if  1 “ dr  G A I 


0 .if  i £ A 
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V 

Thus,  a is  formed  by  replacing  the  l's  in  z(A)  with  the  entries  from 
the  corresponding  row  of  E.  Equivalently,  one  may  consider  the  matrix  , 
E to  be  augmented  by  inserting  zero-columns  where  zeros  occur  in  z(A). 

V 

The  rows  of  the  resulting  H XL  matrix  are  the  sequences  a . The 

correlation  coefficients  for  the  set  A(A;E)  are  easily  obtained  from 

-T 

the  entries  in  E and  EE  . For  T =£  0,  cv^(t)  “ ev^  eXj  ^ i = d^ 

and  J = i + T = dg,  while  cv^(t)  = 0 if  T does  not  occur  as  a 

difference  d - d . For  V =£  A,  one  obtains 

r s 

cvA(0)  ai  *1 
i=0 

K 

"]£  evr  5Xr  • 
r=l 

which  is  the  vAc  entry  in  EE  . 

For  n a positive  integer  let  F^  denote  the  Fourier  matrix  with 

Vi1*1  entry  equal  to  u/V  1 < v,  i < n,  where  u is  a prim- 

itive n**1  root  of  1.  H will  denote  an  n X n Hadamard  matrix  of 

n 

+l's  and  -l's,  having  only  +l's  in  the  first  column.  We  let  F° 

and  H®  be  the  n X (n-1)  matrices  resulting  from  deleting  the  first 

column  from  F and  H , respectively.  From  the  fact  that 
T n n’  T T 

F F ■ 11  H « nl  , it  follows  that  F (F°)  = H°(H°)  = nl  - J , 

nnnnn  nn  nn  nn 

where  I and  J are  the  n X n identity  matrix  and  the  n X n matrix 
n n 

of  +l's,  respectively. 


3.  THE  TERNARY  AND  COMPLEX  SEQUENCES 


For  each  of  our  sets  A * A(A;E),  A will  be  one  of  two  algebraically 
distinct  types  of  (L,K)  partial  difference  sets,  and  E will  be  a^F®+j  , 
vw  or  a^H^,  where  ■ K The  scalar  normalizes  the  rows 


K+r  K K’ 

of  the  resulting  E's. 


In  each  case  every  entry  of  E has  magnitude  K 
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Thus,  the  correlation  coefficients  for  T =/=  0 all  have  magnitude  1/K 

-T 

or  0.  Every  off-diagonal  entry  of  EE  is  0 for  E = a H , and  -1/K 

lv  K. 

for  E = aKF^+i  or  aKHK+l’  Therefore,  | cv^ (0)  | < 1/K  for  V =£  X. 

Whenever  E is  one  of  these  three  matrices,  we  have  c (A)  **  1/K. 

max 

For  E = aKHK+l  °r  aKHK’  resu^tin8  sequences  are  ternary  with  -a^, 

0,  and  +a^  as  entries.  The  sequences  are  complex-valued  when 

E “ °KFK+1 * 


For  a given  A,  setting  E equal  to  a F°  or  .OH®,  yields  one 
more  sequence  than  setting  E equal  to  a^H^,  since  F°+1  and  H£+1  are 
(K  + 1)  X K matrices,  while  is  a K X K matrix.  The  larger  set  can 
always  be  obtained  since  F^  exists  for  all  n.  However,  for  some  appli- 
cations it  may  be  desirable  to  have  ternary  rather  than  complex  sequences. 
This  may  not  be  possible  with  this  method  since  Hadamard  matrices  Hn  do 
not  exist  for  n not  a multiple  of  4,  with  the  exception  of  n = 2. 


The  first  type  of  partial  difference  set  is  the  planar  (L,K)  set 
with  L >=  K -K+l,  K-l  a prime  power.3'1*  For  K = 2,  we  have 

L ■ 3 and  A = {0,l}.  Letting  E = 2 T®  gives  the  three  sequences 

a1  = 2-,S(l,  1,  0) 
a2  - 2_*S0‘>,  uj2,  0) 
a3  «*  2_,S(u)2,  w,  0) 


2 

where  w is  a cube  root  of  1 satisfying  u>  + o>  + 1 = 0.  The  auto- 
correlations and  cross-correlations  for  the  set  of  three  sequences  are 

a1  * a1  - |(2,  1,  1) 

2*2  1,„  2. 

a * a - (a,  ui. ) 

3*3  1/,  2 . 

a * a - “2(2,  a)  , oj) 

1*2  1.  , 2 . 

a * a - jj(-l,  <*>  , u) 

2*3  1,  . 2. 

a * a - 2<“2.  “>  ) 

2 3 1 

a*  * aJ-  |(-1,  1,  1). 


r 
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Therefore,  this  set  achieves  the  bound  W(3,3)  = 1/2.  This  set  of 
sequences  (for  K = 2)  is  the  smallest  one  guaranteed  by  the  following. 

THEOREM  1.  If  K - 1 is  a prime  power,  then  there  exists  a 

set  of  K + 1 distinct  complex-valued  sequences  of  period 

L = - K + 1 with  c equal  to  the  bound  W(L,K+1)  = 1/K. 

max 

The  sequences  of  Theorem  1 are  constructed  from  the  planar  difference 
sets  for  the  prime  power  K - 1,  and  the  truncated  (K  + 1)  X K Fourier 
matrix  aRF»+r 

The  second  planar  set  occurs  for  K = 3,  L = 7,  and  A = {0,1,3}. 

We  may  let  E = 3 ^F°  or  3 , where 

1111 
1 1-1-1 
1-1  1-1 
Ll  -1  -1  1J 


and 


H? 


“1  1 1" 

1 -1  -1 
-1  1 -1 
bl  -1  U . 


For 

— *5 

E *=  3 H the  four  sequences 

are 

a1  - 3~*(  1,  1, 

0,  1, 

0, 

o. 

0) 

a2  - 3~'*i  1,-1, 

0 

1 

0, 

o, 

0) 

a3  - 3‘Vl,  1, 

0,-1, 

o, 

o, 

0) 

a4  - 3~\-l,-l, 

o,  1, 

o, 

o. 

0). 

The 

bound  W(7,4)  - 1/3 

is  met  by  this 

// 


1 
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The  only  known  planar  difference  sets  are  those  for  K - 1 a prime 
power.  Since  2 is  the  only  prime  power  congruent  to  2 modulo  4,  the 
four  sequences  above  form  the  only  ternary  set  resulting  from  this  con- 
struction which  achieve  the  bound. 

For  K - 1 = 3,  the  set  A = {0,1, 3, 9}  is  a planar  difference  set 

modulo  13.  For  this  set,  the  matrix  ■=■  H,  can  be  used  to  construct 

2 4 

four  sequences 

a1  - ^(1,  1,  0,  1,  0,  0,  0,  0,  0,  1,0,  0,  0) 

a2  = |(1,  1,  0,-1,  0,  0,  0,  0,  0,-1,  0,  0,  0) 

a3  «=  i(l,-l,  0,  1,  0,  0,  0,  0,  0,-1,  0,  0,  0) 

aA  *=  |(1,-1,  0,-1,  0,  0,  0,  0,  0,  1,  0,  0,  0). 

The  correlation  coefficients  (excluding  the  auto-correlation  peaks)  are 

all  - 0,  or  i.  The  Welch  bound  for  this  case  is 

1 

W(13,4)  = (— ) = 0.243.  This  set  is  the  smallest  of  the  type  guaranteed 

by 

THEOREM  2.  If  K - 1 is  a prime  power  congruent  to  3 modulo  4y 
then  there  exists  a set  of  K distinct  ternary  sequences  of  period 
L ■=  lP  - K + 1 which  nearly  achieve  the  bound.  In  particular 
= 1/K,  while  W(L,K)  = + l)'*. 

max 

The  sets  of  Theorem  2 use  E = Hadamard  matrix 

Theorem  2 may  be  constructed  from  the  quadratic  residues  in  the  finite 
field  GFCK-l).3*4 

For  the  second  type  of  partial  difference  set  K may  be  any  prime 
2 

power,  and  L = K - 1.  These  are  special  cases  of  the  relative  differ- 
ence sets  of  Elliott  and  Butson.5  A maximal  linearly  recurring  sequence 

2 

of  degree  2 over  the  finite  field  GF(K)  has  period  L = K - 1.  Every 
non-zero  member  of  CF(K)  occurs  exactly  K times  in  such  a sequence, 
while  0 occurs  K - 1 times.  The  set  A of  positions  at  which  the  l's 


»K  f°r 
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occur  form  an  (L,K)  partial  difference  set.  For  K “ 2 this  set  corres- 
ponds to  the  planar  set  in  Z„.  For  K = 3,  the  linearly  recurring 

J 2 

sequence  associated  with  the  primitive  quadratic  x + x + 2 is 
(1, 1,0, 1,2, 2, 0,2) . The  resulting  relative  difference  set  is  A = {0,1,3} 
in  Zg.  The  following  two  theorems  result  from  the  relative  difference 
sets  with  K a prime  power. 

THEOREM  3.  If  K is  a prime  power y then  there  exists  a set  of 

K + 1 distinct  complex-valued  sequences  of  period  L = - 1 

with  c = 1/K.  If  K = 3 (mod  4),  then  a set  of  ternary 
max 

sequences  with  the  same  parameters  exists. 

THEOREM  4.  If  K is  a power  of  2}  then  there  exists  a set  of 

K distinct  ternary  sequences  of  period  L = - 1 with 

c = 1/K. 
max 

The  ternary  sequences  of  the  second  part  of  Theorem  3 are  constructed 
using  the  quadratic  residue  Hadamard  matrix  over  GF(K),  that  is, 

E = Cij,H°+i . The  Hadamard  matrix  H^,  required  in  Theorem  4 can  be  con- 
structed by  the  iterated  tensor  product 


n 


Li  -u 


H„  = H_  ® H , m = 2r  > 2,  r = 1,2,3...  . 
zm  z m 

The  sets  of  Theorems  3 and  4 approach  the  bound  W(L,K)  as  K and  L 
increase.  The  bounds  are 

W (L , K+l ) = (K/(K3  +.K2  - K - 2))**  for  Theorem  3, 


W(L,K)  = ( (K-l)  / (K3  -.K  - l))*5 


for  Theorem  4. 


Clearly  each  of  these  bounds  approaches  1/K  es  K increases.  Thus,  the 
sets  for  Theorems  2,  3,  and  4 are  "asymptotically  optimal." 
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4.  THE  BINARY  SEQUENCES 

In  this  section  we  construct  sets  of  sequences  of  period  L = 2^m  - 1 
from  linear  cyclic  error-correcting  codes  over  GF(2)  of  length  L. 

For  y an  element  of  GF(2n)  let  u(Y)  denote  a linearly  recurring 
sequence  associated  with  the  minimal  polynomial  of  y over  GF(2).  The 
period  tt(y)  of  u(y)  is  the  multiplicative  order  of  y in  GF(2n),  and 
must  divide  2n  - 1.  Let  v(y)  be  the  (2n  - 1) -dimensional  vector  over 
GF(2)  formed  by  juxtaposing  l(y)  periods  of  the  sequence  u(y),  where 
I(Y)  = (2n  — 1) / 7i (y)  • If  Y is  a primitive  element  of  GF(2n) , then 
I(Y)  = 1>  and  v(y)  together  with  all  of  its  cyclic  shifts  generate  an 
n-dimensional  cyclic  code  over  GF(2)  with  minimum  weight  2°  ^ . Now 
suppose  6 = yr»  and  let  C(r)  be  the  code  generated  by  all  cyclic  shifts 
of  both  v(y)  and  v(6).  C(r)  is  a cyclic  (n+d) -dimensional  code,  if 
6 is  not.  conjugate  to  y,  where  d is  the  degree  of  6 over  GF(2).  Kasami, 
Lin,  and  Peterson  have  computed  the  weight  distributions  of  several  such 
codes.6  C(r)  decomposes  into  classes  of  cyclically  equivalent  vectors. 

By  selecting  one  representative  from  each  class  of  period  2n  - 1 and 
replacing  each  0 and  each  1 with  -(2°  - 1)  **  and  +(2n  - 1) 
respectively,  a set  A(r)  of  normalized  binary  sequences  is  obtained. 

The  Gold  sequences  result  from  letting 

2(n+l)/2  + x for  n Qdd 

r “ 

,(n+2)/2  . . . 

2 +1  for  n even 

In  this  case  C(r)  has  dimension  2n,  and  the  2^n  - 1 non-zero  code 
words  fall  into  2n  + 1 classes,  each  of  period  2n  - 1,  provided  r is 
relatively  prime  to  2n  - 1-.  (When  g.c.d.  (r,  2°  - 1)  «*  s > 1,  C(r) 

contains  only  2n  classes  of  period  2n  - 1.  The  remaining  2n  - 1 non- 
zero code  words  have  period  (2n  - l)/s.)  The  family  A(r)  contains 

H ■ 2°  + 1 sequences  of  period  L » 2n  - 1,  with  c ■ r/L  (see 

max  , 

footnote  2).  Thus,  for  the  Gold  sequences  cmax  is  approximately  2L  , 

2**L  **  for  n even,  odd,  respectively;  whereas  W(L,M)  «■  W(2n  - l,2n  + 1) 
_L 

is  approximately  L . 
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For  our  sequences  we  let  n = 2m  and  r = 2®  + 1.  The  resulting 

code  C(r)  has  length  L = 2^®  - 1,  dimension  3m  and  minimum  weight 

2m  *(2®  - 1).  In  fact  the  only  non-zero  weights  occurring  in  C(r)  are 

2®  3 (2®  - 1),  2^®  ^ , and  2m  3 (2®  + 1),  see  footnote  6.  Thus,  the 

cross-correlation  coefficients  and  auto-correlation  side-lobes  for  the 

set  A(r)  of  sequences  of  period  L assume  only  the  three  values 

(-2®  - 1)/L,  -1/L,  and  (2m  - 1)/L.  It  follows  that  c = (2®  + 1)/L 

m max 

“ 1/(2  - 1).  This  approximates  c for  Gold  sequences  of  period 

2ro+l  max 

2 ~ 1 “ 2L  + 1.  The  3m-dimensional  code  C(r)  contains  a single 

class  of  vectors  of  period  2m  - 1 generated  by  6.  The  remaining 

23m  _ 2m  vectors  fall  into  2®  classes  of  period  L = 2^®  - 1.  Thus, 

M = 2®  and  the  bound  is  W(L,M)  = ((2®  - l)/(23®  - 2®  - I))*5!  which  is 

-Jj 

approximately  equal  to  L 


The  first  example  of  this  construction  is  for  m = 2,  L = 15,  M = 4. 

4 5 

Letting  y be  a root  of  x + x + 1,  and  6 = y , we  have 

v(Y)  =(000100110101111), 

v(6)  =(011011011011011). 

The  four  cyclically  distinct  vectors  in  C(r)  are  v(y)  together  with 

v(y)  added  to  each  of  the  three  shifts  of  v(6) . For  the  associated 

set  A(r)  we  have  L = 15,  M = 4,  and  c = 1/3;  while 
« max 

W(15,4)  = (3/59K  = 0.225. 
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